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Abstract 



We review the composition of Jupiter-family comet dust as inferred from infrared spectroscopy. We find that Jupiter-family 
\l comets have 10 (im silicate emission features with fluxes roughly 20-25% over the dust continuum (emission strength 1.20-1.25), 
similar to the weakest silicate features in Oort Cloud comets. We discuss the grain properties that change the silicate emission 
feature strength (composition, size, and structure/shape), and emphasize that thermal emission from the comet nucleus can have 
("^ significant influence on the derived silicate emission strength. Recent evidence suggests that porosity is the dominant parameter, 
O ,3.1 though more observations and models of silicates in Jupiter-family comets are needed to determine if a consistent set of grain 
' parameters can explain their weak silicate emission features. Models of 8 m telescope and Spitzer Space Telescope observations 
^~ have shown that Jupiter-family comets have crystalline silicates with abundances similar to or less than those found in Oort Cloud 
comets, although the crystalline silicate mineralogy of comets 9P/Tempel and C/1995 Ol (Hale-Bopp) differ from each other in Mg 
and Fe content. The heterogeneity of comet nuclei can also be assessed with mid-infrared spectroscopy, and we review the evidence 
i i for heterogeneous dust properties in the nucleus of comet 9P/Tempel. Models of dust formation, mixing in the solar nebula, and 
comet formation must be able to explain the observed range of Mg and Fe content and the heterogeneity of comet 9P/Tempel, 
' ^ although more work is needed in order to understand to what extent do comets 9P/Tempel and Hale-Bopp represent comets as a 
JJT whole. 
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1. Introduction 

The dust in Jupiter-family comets (JFCs) is considered 
to be pristine, that is, little or no processing of the dust has 
occurred since incorporation into the comet nucleus. We 
do not equate pristine dust with pre-solar dust. The pres- 
ence of crystalline silicates in the dust comae of JFCs and 
Oort Cloud (OC) comets, and the absence of crystalline 
silicates in the interstellar medium (< 2.2%; Kemper et 
al., 2004, 2005), strongly suggests that much of the dust 
in comets has been thermally processed (T gl 1000 K). 
Strong radiogenic heating may occur in comets of suffi- 
ciently large (radius, R n , > 10-20 km) porous and dust rich 
nuclei (Prialnik and Podolak, 1999; Prialnik ct al., 2008; 
Merk and Prialnik, 2006; McKinnon et al., 2008), but the 
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presence of highly volatile molecules, such as CO and/or 
CO2 observed in the comae of 9P/Tempel (A'Hearn et al., 
2005), 21P/Giacobini-Zinner (Mumma et al., 2000), and 
103P/Hartlcy (Crovisier et al., 2000), suggests that comet 
interiors formed at very low temperatures (as low as « 30 K 
for CO to condense) and have since remained cold (Ya- 
mamoto, 1985; Crovisier, 2007). The simultaneous pres- 
ence of low-temperature and high-temperature solar neb- 
ula condensates in comet nuclei is possible if thermally pro- 
cessed dust can be mixed from the hot inner-solar nebula 
into the cold comet forming zone (heliocentric distance, r^, 
^ 4 AU) perhaps by turbulent diffusion (Wehrstedt, 2003; 
Wehrstedt and Gail, 2003, 2008; Wooden ct al., 2007) or by 
large-scale radial (meridional) flows (Keller and Gail, 2004; 
Ciesla, 2007; Tscharnutcr and Gail, 2007). 

Mid-infrared (mid-IR) spectroscopy (A ~ 10 /im) is an 
effective technique for assessing the composition of comet 
dust in a large number of comets: 1) the peak thermal 
emission output from comet dust at r - /, w 0.5-4.0 AU is 
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at A w 7 — 20 /im; 2) the mid-IR covers prominent emis- 
sion features from silicate minerals, a major constituent of 
comet dust; and 3) only low to moderate spectral resolu- 
tion (A/AA > 100) is required to critically assess the com- 
position of comet dust. By studying the physical properties 
of dust ejected by comet nuclei, we gain insight into the 
formation mechanism of dust in the solar nebula, and the 
structure and dynamics of the proto-planetary disk (PPD). 

Through comparisons between dust in the interstel- 
lar medium, interplanetary dust grains (Brownlee, 1985), 
comet dust, and asteroidal materials, we may attempt 
to ascertain the extent to which PPD processes altered 
pre-solar dust and mixed it throughout the comet forming 
zone (Wooden, 2008). Being a mixture of pre-solar and so- 
lar nebula processed material, measuring the composition 
of comet dust helps us understand both the origin of the 
constituents of the primitive solar system, and their sub- 
sequent evolution into planetcsimals (Ehrcnfrcund ct al., 
2004; Wooden, 2008). In a similar vein, we can use JFCs 
to learn about the compositions of Scattered Disk objects 
since these bodies are dynamically linked (Duncan et al., 
2004) and therefore their dust properties may be related. In 
contrast, Kuiper Belt objects from the classical disk likely 
formed in situ, and represent a population of icy bodies 
separate from the Scattered Disk objects and comet nuclei. 
Indeed, Kuiper Belt objects have redder surface colors, on 
average, than comet nuclei (Jewitt, 2002), indicating that 
they have a different origin or evolutionary history. 

Since mid-infrared spectroscopy is effective at measuring 
comet dust properties, it offers important insight into the 
structure of the PPD. We can use mid-IR observations of 
comet dust to explore the heterogeneity of the PPD on sub- 
kilometer and super-kilometer (i.e., possibly tens of AU) 
scales by measuring the heterogeneity of comet nuclei. Dust 
may be ejected into the comet coma from a localized re- 
gion on the nucleus at different epochs either by natural 
or unnatural processes, e.g., seasonal variations in insola- 
tion, outbursts, or collisions (Deep Impact-like events or 
natural). If the dust composition of the ambient coma dif- 
fers from the localized regions, then we may conclude that 
the comet nucleus is heterogeneous. Nucleus heterogene- 
ity could arise from the aggregation of varied planetesimals 
(see, e.g., Belton et al., 2007) and evidence for heterogene- 
ity implies the proto-planetary disk in the comet forming 
regions was heterogeneous on kilometer to sub-kilometer 
sized scales (the typical radius of a comet nucleus is of or- 
der 1 km; Lamy et al. , 2004) . Alternatively, the dust prop- 
erties of individual nuclei could be relatively homogeneous, 
yet the dust compositions vary from comet-to-comct. This 
second scenario favors PPD heterogeneities on the super- 
kilometer scale. A taxonomy of comets classified according 
to their dust properties has the promise to provide evidence 
supporting local (spatial or temporal) variations in dust 
composition in the comet forming region of the PPD. Ev- 
idence for chemical heterogeneity on both sub- and super- 
kilometer scales has been discovered through observations 
of comet gas species (Bockelee-Morvan et al., 2004), e.g., 



the C2 comet classes (A'Hearn et al., 1995), the organic 
chemistry classes (Mumma et al., 2003; Bockelee-Morvan 
et al., 2004; Crovisier, 2007), the apparent CO variability of 
comet 21P/Giacobini-Zinner (Mumma et al., 2000; Weaver 
et al., 1999), and the CO2 and H 2 composition of comet 
9P/Tcmpcl (Feaga et al., 2007). 

In this paper, we review the composition of Jupiter- 
family comets as inferred from infrared spectroscopy: in §2 
we review mid-IR spectroscopic observations of JFCs; in §3 
we summarize and discuss observations of silicate mineral- 
ogy in JFCs, we demonstrate (§3.3) that the silicate feature 
strength is increased by subtracting emission from the nu- 
cleus, but not increased enough to match the strong silicate 
features of some OC comets, and from this perspective we 
discuss opportunities for important observations in the fu- 
ture; in §4 we review the mid-IR observations and models 
of comet 9P/Tcmpel in the context of the Deep Impact mis- 
sion; in §5 we discuss the potential significance to our un- 
derstanding of PPD processes of the detection of hydrated 
and Fe-rich minerals identified in comet 9P/Tcmpcl; and, 
in §6 we summarize the review. 

2. Mid- infrared Spectroscopic Observations 

Mid-infrared spectra, with A/AA ^ 50, exist for 13 
Jupiter-family comets and these observations are summa- 
rized in Table 1 . We also list whether or not a silicate emis- 
sion feature was detected above the continuum for each 
spectrum, and the strength of that feature, if available. It 
is important to note that the 10 /jm continuum in a mid-IR 
spectrum of a comet is composed of thermal emission from 
dust species that exhibit little variation in emissivity with 
wavelength for A = 7 — 14 /im. The common "continuum" 
contributers are amorphous carbon, FeS, and grains large 
enough to be optically thick to 10 /im emission (Hanner 
et al., 1994; Wooden, 2002). Thermal emission from the 
nucleus can also be an important source of continuum 
flux, and should be removed before measuring the silicate 
emission (see §3.3). A silicate feature has been detected in 
8 of the 13 comets in Table 1, and the typical silicate emis- 
sion band strength is 1.20-1.25. The silicate emission band 
strength is computed from the ratio of the mean in-band 
flux between 10.0 and 11.0 ^m to a scaled Planck function 
that is fit to small sections of the smooth continuum at 
A < 8.2 - 8.4 11m and A > 12.4 - 12.5 ^m and evaluated 
at 10.5 ptm (Hanner et al, 1996; Sitko et al., 2004; Sugita 
et al., 2005). Observed silicate band emission strengths in 
Oort Cloud comets range from 1.0 to 4.0 (Hanner ct al., 
1994; Williams et al., 1997; Harkcr et al., 1999; Sitko et 
al. , 2004) . Such strong silicate band emission has not been 
observed in a JFC, except in the post-Deep Impact coma 
of 9P/Tcmpcl (Lisse et al. 2006; Harker et al. 2007; §4). 

The Spitzer Space Telescope (Werner et al., 2004; Gehrz 
et al., 2007) has significantly increased the number of spec- 
troscopic observations of JFCs in the mid-IR. Three comets 
with high-quality spectra are listed in Table 1, and 9 more 
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are briefly summarized by Kelley (2006). This boon to 
JFC spectroscopy is largely attributable to the high sen- 
sitivities of the Infrared Spectrograph (IRS) instrument 
(Houck et al. , 2004) and the high observing efficiency of the 
Spitzer spacecraft. As future infrared space telescopes will 
have larger primary mirror diameters, such as the JAXA 
Space Infrared Telescope for Cosmology and Astrophysics 
(SPICA), and the NASA James Webb Space Telescope, 
bright comets may become unaccessible to space-borne in- 
strumentation. Instead, ground-based and airborne tele- 
scopes will be used to observe the brightest comets. 

Ground-based observations account for 1/2 of the spec- 
tra in Table 1; about 2/3 of those spectra were taken 
at the NASA's 3 m Infrared Telescope Facility (IRTF). 
Eight meter telescopes hold great promise for the spectro- 
scopic investigation of JFCs. The Gemini and Subaru tele- 
scopes have already provided three definitive detections of 
crystalline silicates in JFCs; the measured crystalline-to- 
amorphous sub-/im mass ratios are: ^0.9 for 78P/Gchrcls 
(Watanabe et al. 2005; also discussed in Ootsubo et al. 
2007a), 0.3-0.8 for 9P/Tcmpcl (Harker et al., 2005, 2007; 
Sugita et al., 2005; Lisse et al., 2006), and qualitatively 
described as less than in Hale-Bopp for 73P/Schwassmann- 
Wachmann (Harker et al., 2006b; Sitko et al., 2006d). Also, 
Gcmini-N and Subaru observations of the Deep Impact 
event provided, in part, the evidence for the heterogeneous 
structure and composition of the nucleus of 9P/Tempel 
(Harker et al., 2007; Kadono et al., 2007). With the con- 
tinued use of mid-IR instruments on large telescopes, and 
in conjunction with the archive of space-based observa- 
tions, we can increase the number of high-quality spectra 
of JFCs to statistically significant numbers, which will give 
us the opportunity to establish the extent of crystalline 
silicate enrichment in JFCs and the degree of JFC nucleus 
heterogeneity. 

Future airborne instrumentation on the Stratospheric 
Observatory for Infrared Astronomy (SOFIA) will cover 
mid-IR wavelengths in spectroscopy mode (Becklin et al., 
2007). With the telescope executing science operations at 
altitudes ft 12.5 km, SOFIA will be an important resource 
for comet astronomy at 15-40 fim, where atmospheric ab- 
sorption greatly interferes with observations at ground- 
based observatories. At these wavelengths, the spectral res- 
onances of crystalline silicates are more separated in wave- 
length, and therefore more distinct, than the resonances 
near 10 ^m. With 15-40 /im spectra, we can accurately 
measure the mincralogical makeup of silicate crystals in 
a comet's coma by measuring the central positions and 
strengths of the crystalline emission resonances. In partic- 
ular, the Fc to Mg ratio of silicate crystals is better con- 
strained when the longer wavelength peaks are measured. 
SOFIA will add vital observations to the few crystalline 
measurements made at these wavelengths. 



3. Silicates in Jupiter-Family Comets 

3.1. Properties of silicate grains and their thermal emission 

Silicate minerals in comets have prominent emission 
features in the mid-IR. Molecular bonds between Si and 
O atoms have stretching and bending modes that pro- 
duce emission at 8-12 /zm and at 15-40 /im, respectively. 
The central wavelengths, relative strengths, and shapes of 
these features are diagnostic of the mineral composition. 
Wooden (2002) reviews the silicate minerals observed in 
ground-based mid-IR spectra of comets: [Mg y ,Fei_ y ] 2 Si04 
(olivine) and [Mg^Fei-JSiOs (pyroxene), where x and y 
range from (Fe-pure: fayalite and ferrosilite) to 1 (Mg- 
pure: forsterite and enstatite). Cometary silicates are found 
in both amorphous (non-stoichiometric) and crystalline 
structures. An example of the narrow thermal emission 
resonances attributed to silicate crystals, and the broad 
emission attributed to amorphous grains is presented in 
Fig. 1. 

Comet comae are optically thin to mid-IR radiation. In 
order to detect the thermal emission from silicate grains in 
a dust coma, the grains must be in sufficient quantities to be 
detected above the observational uncertainty of the smooth 
thermal continuum (Hanner et al., 1994). Aside from the 
amount, each grain's size, composition, and structure play 
important roles in the strength of the silicate emission, pri- 
marily by determining the grain equilibrium temperature. 
We discuss each of these properties (size, composition, and 
structure) in turn. 

First, dust grain size affects the grain temperature. Sub- 
orn sized grains have higher temperatures than (solid) 
grains larger than 1 ,um of the same composition. The 
temperature difference arises from the dependence of grain 
emission efficiency on grain size. All silicate dust grains 
greater than 0.1 /zm in size are efficient absorbers of sun- 
light at visible wavelengths (the peak of the solar output), 
but the sub-/xm sized grains are inefficient emitters in the 
mid-IR. Since grain radiative cooling at r] x m 0.5 — 4 AU is 
dominated by mid-IR radiation, the temperatures of sub- 
orn sized grains are warmer, which subsequently increases 
their total thermal emission flux. 

In contrast with grain size, which primarily affects the 
emission efficiencies of comet grains, composition affects 
grain absorption efficiencies. Increasing the amount of Fe 
in a silicate grain (decreasing Mg) increases the grain's 
absorption efficiencies at visible and near-IR wavelengths, 
which in turn increases the grain's equilibrium temperature 
(Dorschner et al., 1995). The higher temperatures cause Fe- 
rich grains to be easier to detect by spectroscopic means, 
yet Harker (1999) found the grain temperatures of amor- 
phous silicates in comet Hale-Bopp were best-modeled with 
roughly equal amounts of Mg and Fe (x = y = 0.5). Fur- 
thermore, Harker (1999) found the silicate crystals to be 
distinctly Mg-rich (x = y = 0.9 — 1.0). These grain pa- 
rameters have been commonly used by other investigations 



3 



(Moreno ct al., 2003; Honda et al., 2004; Min et al., 2005; 
Harkcr ct al. , 2007) . However, Lisse et al. (2006) found sig- 
nificant amounts of Fe-rich pyroxene in the Deep Impact 
ejecta from comet 9P/Tempel (see §4 and Fig. 4). In Fig. 2, 
we summarize the variation of the peak wavelength position 
with Mg content for strong crystalline olivine and ortho- 
pyroxene peaks observed near 10 /xm from the analyses of 
Koike et al. (2003) and Chihara et al. (2001). 

Finally, grain structure (i.e., porosity and shape) influ- 
ences the thermal emission from dust. For two otherwise 
similar grains, increased porosity will produce a warmer 
equilibrium temperature and, in larger grains, a stronger 
10 fim silicate emission feature (Hage and Greenberg, 1990; 
Xing and Hanner, 1997; Min et al., 2006; Kolokolova et 
al., 2007; Voshchinnikov and Henning, 2008). The effect 
is due to insulating vacuum inside porous grains, caus- 
ing sub-domains to be increasingly thermally isolated from 
each other. If the porosity of an aggregate grain is more 
than 63% then the temperature begins to approach that of 
the monomer (^0.1 ^m sub-grains): at 80% porosity the 
temperature of the aggregate is only 4-6% lower than a 
monomer (Xing and Hanner, 1997). Also, the shape of a 
porous aggregate determines whether larger porous grains 
lose or retain their silicate resonant features (Kolokolova 
et al., 2007). The net result of increased porosity is a ther- 
mal emission spectrum with a hotter temperature and a 
stronger silicate feature compared to a similar sized com- 
pact grain. Thermal emission models for porous aggre- 
gates that include crystals (computed with the Discrete 
Dipolc Approximation; Moreno et al., 2003) yield similar 
crystalline fractions as the less time-intensive approach of 
modeling a grain distribution composed of discrete crystals 
and mono-mineralogic porous grains (Harker et al., 2002, 
2004a). 

The above three grain parameters (size, composition, and 
structure) all affect the observed strength of a coma's sili- 
cate emission feature, but each parameter also affects the 
shape of the feature. When we consider the silicate emis- 
sion feature's shape, and the shape of the entire spectrum, 
we find that there are no one-to-one degeneracies among 
the dust grain parameters, although there are consistent 
trends. For example, the silicate feature strength can be 
raised by increasing the porosity of grains above ~ 80% vac- 
uum, which yields silicate features and temperatures simi- 
lar to their monomers (Xing and Hanner, 1997). Increasing 
the porosity of dust grains in a coma also shifts the spec- 
tral slope of the Rayleigh-Jeans side of the observed spec- 
trum towards the blue (i.e., there are fewer, cool grains con- 
tributing to the flux at long wavelengths). A small portion 
of this change in spectral slope can be accounted for by flat- 
tening the size distribution to increase the relative number 
of larger, cooler grains. The entire effect, however, cannot 
be accounted for in this manner — this is evident when we 
consider that the largest, highly porous grains have higher 
equilibrium temperatures than the largest grains with more 
compact shapes. Oort cloud comet Hale-Bopp is an exam- 
ple of a comet that displays an extremely strong silicate 



feature, and that fitting its spectrum requires a size distri- 
bution that possesses both highly porous grains and very 
small grains (peak grain radius of 0.15-0.2 /im; Harker et 
al., 2002; Moreno ct al., 2003; Min ct al., 2005). Wooden 
(2002) discusses how comet Halc-Bopp's grain porosity and 
grain size distribution slope are each constrained by spec- 
tra. In another example, Oort Cloud comet C/2001 Q4 
(NEAT) has a silicate feature identical in shape to, but 
weaker in strength than, that of Hale-Bopp (Wooden et al., 
2004). In fact, Q4's silicate feature strength varies on short 
times scales and is best explained by a variable silicate-to- 
amorphous carbon mass ratio (Wooden et al., 2004; Harkcr 
et al., 2004b). Specifically, decreasing the silicate feature 
strength by a factor of 1.5 by removing the grains smaller 
than 0.6 /im broadens the shape of the silicate feature more 
than what is observed for Q4. Therefore grain size effects 
cannot explain the observed variations in the spectra of 
comet Q4. 

With the above arguments in mind, we can consider the 
positive detection of silicate features in many JFCs (Ta- 
ble 1) and draw some conclusions. The weak silicate band 
emission in JFCs, as compared to some OC comets, could 
be due to: 1) silicates may be less abundant with respect 
to other dust components (e.g., amorphous carbon, FeS) in 
JFCs than in OC comets; 2) silicate grains may be larger in 
JFCs than in OC comets; 3) JFC dust may be less porous 
than OC comet dust; or 4) a combination of 1, 2, and/or 
3. Out of the few JFCs that have been studied in-detail, 
there is no clear agreement in the silicate dust content of 
their comae. Kelley et al. (2006) found Jupiter-family comet 
2P/Encke had a peak grain size distribution of 0.4-0.5 /im 
and a sub-/xm silicate dust mass fraction of < 31% (3a 
upper-limit). Hanner et al. (1996) reasonably reproduced 
the silicate feature of 19P/Borrclly to within 10% using a 
50-50 mixture of silicates and carbon. Models of prc-Deep 
Impact spectra of 9P/Tempel favored a silicate dust dis- 
tribution with a peak grain size of 0.7-0.8 /im for solid to 
somewhat porous grains (Harker et al., 2007). It may be 
that these three comets are truly different in their composi- 
tions and grain size distributions, however the models could 
be refined with new or improved constraints. For example, 
the comet Encke model fits could benefit from higher sen- 
sitivity spectra at A < 14 /im that better constrain the 
grain temperatures and silicate compositions, and all three 
comet models could be improved with constraints derived 
from scattered light observations (§3.2). Thermal emission 
modeling has the potential to constrain grain parameters, 
but only within the limitations of a given data set. 

A few OC comet thermal emission spectra have also been 
modeled. Models of comets C/2001 Q4 (NEAT), C/2002 V2 
(NEAT), and Hale-Bopp reveal comae with smaller peak 
grain sizes (0.1-0.5 /j,m) that are dominated by crystalline 
silicate grains (see Ootsubo et al., 2007a, and references 
therein). Hanner et al. (1994) find that weak silicate emis- 
sion in several OC comets is correlated with a weak dust 
scattering continuum and weak 3 ^m thermal emission, al- 
together suggesting a lack of small grains. From the thermal 
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emission studies it is not clear which of the possible causes 
(1, 2, or 3) manifest the differences between comet mid-IR 
spectra, but, from the limited results above, OC comets 
seem to have an abundance of the smallest grain sizes, when 
compared with JFCs. In the next section we present recent 
developments in models of light scattered by comet dust 
that suggest grain porosity is the dominant property. 

3.2. Grain Porosity, Thermal Emission, and Light 
Scattered by Comet Dust 

Models and observations of light scattered by comet dust, 
in tandem with mid-IR observations, are providing intrigu- 
ing evidence that supports grain structure as the domi- 
nant cause for the difference between JFCs and OC comets. 
In light scattering models, grain structure refers to one of 
three types: solid grains; compact porous aggregates, de- 
scribed by ballistic particle-cluster aggregation (BPCA); 
or, fluffy (highly porous) aggregates, described by ballistic 
cluster-cluster aggregation (BCCA) (Mukai et al., 1992). 
Grain structure affects the strength of the silicate feature, 
polarization, and surface brightness radial profiles (the lat- 
ter due to the changing efficiencies of gas drag and the so- 
lar radiation force; Meakin and Donn, 1988; Mukai et al., 
1992; Nakamura et al., 1994). 

Kolokolova et al. (2007) delineate two types of dust co- 
mae: 1) comae dominated by compact grains that are either 
solid, or compact porous; or 2) comae dominated by fluffy 
aggregates. In scattered light models (Kimura et al., 2003; 
Kimura and Mann, 2004; Kimura et al., 2006), both grain 
types reasonably reproduce the observed optical properties 
of comet dust (for a review, see Kolokolova et al., 2004), 
but one observable breaks the degeneracy: the relationship 
between the 10 /xm silicate emission and the structure of 
the aggregates. 

Kolokolova et al. (2007) show that the surface brightness 
distributions of the dust and gas in comet comae vary from 
comet-to-comet. They suggest that gas-rich comets have 
dust comae that are more spatially compact than their gas 
comae (e.g., comet Encke, Jockers et al., 2005), and dusty 
comets have dust and gas comae with similar surface bright- 
ness distributions. Gas-rich comets have lower optical po- 
larization maxima than dusty comets, and the low optical 
polarization is correlated with weak or absent 10 [im sili- 
cate features (Levasseur-Rcgourd et al., 1996). Exploiting 
the correlation between low optical polarization and dust 
surface brightness distribution, and using recent observa- 
tions that have shown the low optical polarization is due 
to gas contamination in broadband filters (Kiselev et al., 
2004; Jcwitt, 2004; Jockers et al., 2005), Kolokolova et al. 
(2007) develop a physical explanation to describe the cor- 
relation between dust-to-gas ratios and the 10 pm silicate 
emission feature. Their argument is as follows: 1) fluffy, 
porous aggregates have stronger or equivalent 10 /*m sili- 
cate emission features compared to compact, or solid, ag- 
gregates (summarized above); 2) fluffy aggregates should 



be more easily accelerated to the gas outflow velocity than 
compact grains, primarily due to their larger surface areas 
per mass (Meakin and Donn, 1988; Nakamura et al., 1994); 
3) in constrast, compact grains are not accelerated to the 
gas expansion velocity, and dust comae consisting of com- 
pact grains will be more concentrated around the comet 
nucleus than the gas. Altogether, the scattered light and 
the thermal emission properties of comet dust suggest that 
grain porosity could be the predominant characteristic that 
determines the strength of the comet's silicate features. 

3.3. The influence of the nucleus on observed silicate band 
strength 

When assessing the dust content of comet comae, it is im- 
portant to remove the contribution of the nucleus from the 
total thermal emission flux. The nucleus can be a dominant 
contributor to the mid-IR flux of JFCs, which are typically 
less active than their Oort Cloud counterparts (A'Hearn 
et al., 1995). A large nucleus flux skews the silicate band 
strength towards lower values. To illustrate this fact, we 
estimate the fraction of the mid-IR flux attributable to the 
nucleus in a typical JFC observed at r] x = 1.8 AU. We 
will base our model comet on surveys of JFC dust produc- 
tion rates (A'Hearn et al., 1995) and nucleus sizes (Lamy 
et al., 2004). We will also "observe" our model comet at 
r/j = 1.0 AU by appropriately increasing the dust produc- 
tion rate and phase angle to show how the nucleus fraction 
varies with heliocentric distance. We then examine Table 1 
to determine if the weak silicate features are actually a 
product of large nuclei with low dust production rates. 

Since visible light observations of comets arc more readily 
available than mid-IR observations, it is advantageous to 
use the visible observations to estimate the coma brightness 
of a typical JFC. The visible light proxy for dust production 
is the quantity Afp, where A is the Bond albedo of the 
dust, / is the dust areal filling factor in the field-of-view, 
and p is the projected radius of the ficld-of-view, typically 
expressed in units of cm (A'Hearn et al., 1984). The filling 
factor may be expressed as / = Nair^ 1 p~ 2 , where Na is 
the total cross-section of the grains in cm 2 . Technically, 
/ depends on the extinction cross-sections of the grains 
and not on the geometric cross-sections (A'Hearn et al., 
1984). For comae with constant isotropic outflows, Afp is 
independent of aperture size. 

The thermal emission, F\ in units of W cm~ 2 ^m -1 , from 
a collection of isotropically emitting isothermal grains is 

2 

F x = (l-A)nB x (T)f^, (1) 

where A is the mean bolometric albedo of the dust (« 0.32; 
Gchrz and Ney, 1992) , B\(T) is the Planck function in units 
of W cm~ 2 pm^ 1 sr _1 evaluated for the grain temperature 
T and at a wavelength of A, and A is the comet-observer 
distance in cm. It is useful to express Eq. 1 in terms of Afp 
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and 9, where 9 — 206265 p A 1 is the angular size of the 
aperture in units of arcseconds, 



Fx = A( } 206265 ~ A - " 



(2) 



Here, we explicitly write the phase angle dependence of 
the albedo as A(a); the thermal emission term, 1 — A, is 
independent of phase angle. Visual albedo measurements 
of Jupiter-family comets are rare. We found that only one 
comet, 21P/Giacobini-Zinner, has been studied, and it has 
a low albedo: A(64°) = 0.07 - 0.15 (Tclesco et al., 1986), 
A(22°) = 0.11 (Pittichova et al., 2008). In the ncar-IR, 
JFCs have albedos equivalent to or lower than Oort Cloud 
comets (Campins et al., 1982; Hanncr and Newburn, 1989). 
Since comets are typically observed at low-to-moderate 
phase angles, we have chosen A(a < 60°) to be 0.15, which 
is slightly under the apparent mean A(a) of several Oort 
Cloud comets (see Fig. 1 of Kolokolova et al., 2007). The 
mid-IR coma flux is fairly sensitive to A(a), therefore mea- 
sured values should be used when possible. Assuming a 
value of 0. 15 introduces a factor of 2 uncertainty in the mid- 
IR flux when we consider the observed range is w 0.1—0.25 
(excluding comet Halc-Bopp; Kolokolova et al., 2004). We 
take the ensemble average dust temperature to be 10% 
warmer than the blackbody temperature, Tbb, at the ob- 
served rh, 



T w 1.10 Tbb 



306K 



(3) 



which is within the typical range (0-30%) observed in many 
comets (Gehrz and Ney, 1992; Hanner et al., 1996). 

The near-Earth asteroid thermal model (Harris, 1998) 
has proven to be sufficient for modeling the thermal emis- 
sion from comet nuclei, although with some modifications. 
Groussin et al. (2004) suggests an IR-beaming parameter 
of 0.85 for comets, based on the surface models of Lagerros 
(1998). In the formalism of the near-Earth asteroid thermal 
model, the beaming parameter contributes to the observed 
effective temperature of the nucleus. The value of 0.756 has 
been favored for large, main-belt asteroids (Lebofsky et al., 
1986), but equivalent or larger values, up to w 1.0, have 
been experimentally determined with infrared observations 
of comet nuclei (Fernandez et al., 2000; Kelley et al., 2006; 
Fernandez et al., 2006; Groussin et al., 2007). We adopt 
0.85 throughout our paper. We also use a geometric albedo 
of 0.04, a typical value for comet nuclei (Lamy et al., 2004) . 

A'Hearn et al. (1995) surveyed the optical properties of 
85 comets, of which 34 belong to the Jupiter-family class. 
Their median observed Afp value for JFCs is 110 cm at 
a median r/j of 1.8 AU. We derive a median nucleus ra- 
dius of 1.8 km for the 34 JFCs using the radii listed by 
Lamy et al. (2004) and updated values for 81P/Wild and 
9P/Tempel from Brownlee et al. (2004) and A'Hearn et 
al. (2005). Table 2 lists the nucleus-to-coma flux ratios at 
10 /xm for various observing parameters of our median JFC 
(Afp = 110 cm, r h = 1.8 AU, R n = 1.8 km). Of course, 



our fictional "median JFC" parameters are not representa- 
tive of any particular JFC, and that spectra of each comet 
should be considered on a individual basis. Table 2 also lists 
estimates for the median JFC observed at = 1.0 AU by 
scaling the Afp value up by a factor of 1.8 2 3 , as suggested 
by A'Hearn et al. (1995). 

A bright nucleus may affect derived coma parameters, the 
most immediate of which is the silicate band strength. It is 
clear from Table 2 that the nucleus flux can comprise a sig- 
nificant amount of the total flux measured in mid-IR spec- 
tra, as high as 80-90%, especially when observing comets 
with low dust production rates (Afp ~ 100 — 200 cm) using 
the narrow-slits (0.4-1.0" wide) on ground-based instru- 
mentation. In this case, if the silicate band strength appears 
to be between 1.10 and 1.30, then the nucleus subtracted 
(corrected) silicate band strength would lie between 1.50 
and 3.0! With a larger slit, the nucleus could be 10-30% of 
the flux. The corrected silicate band strength would then 
increase to 1.11-1.43. The latter small changes in the sili- 
cate band emission of a JFC may not be significant when 
comparing a JFC to comet Halc-Bopp, which had a silicate 
emission feature that ranged from 3.0 to 4.0, but are im- 
portant when assessing the properties of JFCs as a whole. 

Removing the nucleus is relevant to grain thermal 
emission modeling, even in comets with no discernible 
silicate emission feature. Comet nuclei have temperature 
distributions unlike collections of isothermal grains (the 
dust coma). Comet nuclei have low thermal inertia (e.g., 
Groussin et al., 2007), therefore, the sub-solar point on a 
slowly rotating nucleus will be hottest, and the tempera- 
ture decreases as the sun-zenith angle (#©) increases (in 
the near-Earth asteroid model, T oc cos 1 / 4 9 Q ). In contrast, 
the thermal propagation timescalcs in small grains are fast 
and they quickly equilibrate to a single temperature — a 
temperature that is cooler than the effective temperatures 
of comet nuclei. Thus, the shortest wavelengths of a ther- 
mal emission spectrum will include the highest fraction 
of flux from the nucleus (e.g., see the reduction of the 
2P/Encke spectra in Kelley et al., 2006). If the nucleus 
emission is not removed, a coma modeler would find an 
increased number of the smallest grains (the grain sizes 
that are constrained by the shortest wavelengths) over a 
best-fit model that excludes the nucleus flux. 

We now refer to Table 1 and discuss our estimates of 
how the thermal emission from each comet's nucleus affects 
the observed silicate band strength. Hanner et al. (1996) 
considers the size of the comet nucleus when computing 
the silicate band strength of comet 4P/Faye, so no correc- 
tion is necessary. Based on images of the nucleus of comet 
19P/Borrclly (Buratti et al., 2004), the 8.0 x 3.2 km pro- 
late nucleus contributes 14-36% of the flux in the spec- 
trum of Hanner et al. (1996). The silicate emission fea- 
ture strength is reported by Hanner et al. (1996) to be 
s=a 1.25. Their plot of the data, however, suggests that 
s=a 1.11 is a more appropriate value. Removing the nu- 
cleus flux increases the strength to 1.13-1.19. Stansberry 
et al. (2004) estimate the silicate band strength of comet 
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29P/Schwassmann-Wachmann to be approximately 1.10. 
We use their model nucleus parameters and derive a cor- 
rected band strength of 1.19. Sitko et al. (2004) observed 
69P/Taylor and find the silicate band strength to be 1.23. 
The radius of 69P is estimated to be 2.1 ± 0.6 km by 
Tancredi et al. (2006), which corresponds to 4-19% of the 
thermal emission in the Sitko et al. (2004) mid-IR spec- 
trum. The result is a minor increase in the the silicate band 
strength to 1.24-1.28. The silicate band strength of comet 
103P/Hartley (1.20) measured by Crovisier et al. (2000) 
requires no correction due to the large aperture of the 
ISOPHOT instrument (24" x 24"), and the small nucleus 
size, R n = 0.71±0.13 km (Groussin et al., 2004). We derive 
a nucleus contribution of < 1% for comet 103P. Finally, 
the observation of 103P/Hartley by Lynch et al. (1998) re- 
quires a minor correction, resulting in an increase of the 
measured silicate band strength from 1.15 to 1.22, which 
matches the ISOPHOT derived silicate band strength. 

4. Comet 9P/Tempel and Deep Impact 

Comet 9P was the primary target of the Deep Impact 
mission — a mission designed to impact a comet nucleus to 
determine its underlying strength, structure, and compo- 
sition. In all ways, comet 9P/Tcmpcl is considered to be 
a typical JFC (A'Hearn and Combi, 2007). The Deep Im- 
pact spacecraft delivered a 370 kg impactor into the nu- 
cleus of 9P/Tcmpcl at a relative velocity of 10.3 km s _1 . 
The fly-by spacecraft observed the event from a distance 
of 500 km (A'Hearn et al., 2005). The impact excavated 
<~ 10 5 — 10 6 kg of ice and dust from the surface and inte- 
rior of the comet (Sugita et al., 2005; Harker et al., 2005; 
Kiippcrs et al., 2005; Keller et al., 2007). 

Some of the most detailed mid-IR spectra of any JFC 
were obtained of comet 9P /Tcmpel in support of the Deep 
Impact mission. Three investigations have been presented: 
1) Spitzer observations of the pre- and post-impact coma 
(approximately 30 and 60 min after impact) at 5-40 /im 
(Lisse et al., 2006); 2) Gemini-N observations of the pre- 
and post-impact coma (approximately 7 min cadence for 
3 hr, and additional spectra 24 hr after impact), primar- 
ily at 8-13 fim (Harker et al., 2007); and 3) a Subaru ob- 
servation of the post-impact coma (3.5 hr after impact) at 
8-13 fiia (Ootsubo et al., 2007b). 

4.1. The Pre-Deep Impact Dust Coma 

The Deep Impact spacecraft studied the ambient coma 
and nucleus of comet 9P and provided excellent context 
for observations from other spacecraft and ground-based 
observatories. Farnham et al. (2007) analyzed images of 
9P/Tempcl taken by the Deep Impact spacecraft and found 
that the pre-impact dust coma has three prominent jets, the 
strongest of which originates near comet 9P's south pole. 
Feaga et al. (2007) observed CO2 and H 2 gas in the coma 
and find that CO2 is preferentially ejected from the south- 



ern pole (correlated with the dust) and that H 2 is pref- 
erentially ejected in the sunward direction (on the north- 
ern hemisphere at the time of observation) . The CO2 and 
H2O asymmetries gives strong evidence for either differen- 
tiation of ices in the nucleus, or primordial heterogeneity 
Comet jets will typically be unresolved in mid-IR obser- 
vations, but it would be interesting to investigate whether 
dust from 9P's discrete sources propagate into separate co- 
herent coma or tail structures (Farnham et al., 2007), and 
if the dust properties of these structures could be assessed 
using the high signal-to-noise spectral map of comet 9P 
taken with the Spitzer /IRS instrument. For the purposes 
of this review, we will assume analyses of mid-IR observa- 
tions of the pre- -Deep Impact coma probe a global average 
of the dust properties of this comet. 

Lisse et al. (2006) present a pre-impact Spitzer spectrum 
of 9P/Tcmpcl, but a detailed analysis of this spectrum is 
not given. In Fig. 3 we present the pre-impact spectrum 
taken by Spitzer on 03 July 2005. A model nucleus has 
been removed, and a silicate feature (strength of s=y 20%) is 
clearly observed at 8-12 /im. We estimate the silicate band 
strength to be in the range 1.19-1.25. 

Harker et al. (2007) describe their pre-impact Gemini 
spectra as "essentially smooth and nearly featureless" with 
a silicate band strength close to 1.0 (no silicate emission), 
although an uncertainty is not presented (the nucleus also 
amounts to 80% of their pre-impact flux) . We estimate the 
silicate band strength to be 1.2±0.2 (after nucleus subtrac- 
tion). A spectrum taken in May 2005, about 1.5 months 
prior to the Deep Impact encounter, did show clear silicate 
emission (we estimate the strength to be 1.50 ± 0.18). In 
May 2005 when dust from the southern jet dominates the 
coma (Farnham et al., 2007), the dust mineralogy is dom- 
inated by solid (zero porosity) amorphous pyroxene sili- 
cates with a peak grain size of 0.7 /xm, whereas in July 2005 
pre-impact the ambient coma is comprised of moderately 
porous (46% vacuum for a 10 /jm grain) amorphous car- 
bon and amorphous olivine silicates with a slightly larger 
peak grain size of 0.8 ttm. Note that for a 10 /im radius 
grain, the silicate feature strength increases from negligible 
to s=a 1.2 as the porosity changes from solid (0% vacuum) 
to moderately porous («50% vacuum) (Harker et al., 2002, 
2007). In summary, JFC comae with low dust production 
rates or measured through small slits may at first appear to 
have "featureless" IR spectra (especially at low signal-to- 
noise ratios) , but modeling after nucleus removal can reveal 
micron-sized solid or 10 /im porous silicate grains, as is the 
case for pre-impact comet 9P (Figs. 2 and 5 of Harker et 
al., 2007). 

4.2. The Post-Deep Impact Dust Coma 

Post-impact mid-IR spectra were dramatically different 
from the pre-impact spectra (Fig. 3). Mid-IR spectra of the 
ejecta show the strong resonances of sub-/im sized crys- 
talline silicates of both olivine and pyroxene grains (Lisse et 



7 



al., 2006; Harker et al., 2007; Ootsubo et al., 2007b). This 
discovery is perhaps the most immediate result of the mid- 
IR investigations. A Spitzer /IRS observation of the 10 /xm 
silicate feature (0.64 hr after impact) is presented in Fig. 2. 
We find the post-impact silicate band emission strength to 
be near 2.0 (both the Gemini and Spitzer observations; see 
column 6 in Table 1), the strongest feature ever observed in 
a spectrum of a JFC. The derived crystalline-to-amorphous 
silicate sub-/im mass ratios range from approximately 0.2 
to 4 (Harker et al., 2005, 2007; Lisse et al., 2006; Ootsubo 
et al., 2007b). The values depend on slit width, orientation, 
and time of observation, which could easily explain the wide 
range of values. The pre-impact crystalline-to-amorphous 
silicate ratio has not been significantly constrained (< 1.6, 
la upper-limit Harker et al., 2007). 

Only in two or three other JFCs, 78P/Gchrcls (Watan- 
abe et al., 2005), 73P/Schwassmann-Wachmann (Harker 
et al., 2006b), and, possibly, 103P/Hartlcy (Crovisier et 
al., 2000), have detections of crystalline silicates been 
reported. We may also include the Centaur-like JFC 
29P/Schwassmann-Wachmann (Stansberry et al., 2004). 
For comparison with 9P, Watanabe et al. (2005) found 
a crystalline-to-amorphous silicate ratio of « 0.9 for 
78P/Gchrels. The crystalline-to-amorphous silicate ratios 
of 78P and the Deep Impact ejecta are significantly greater 
than that found in the interstellar medium, < 0.022, (Kem- 
per et al., 2004, 2005), and similar to or smaller than the 
crystalline-to-amorphous ratios measured in the comae 
of some Oort Cloud comets (2-3 for comets Hale-Bopp, 
C/2001 Q4 (NEAT), and C/2002 VI (NEAT); Harker 
et al., 2002, 2004a; Wooden et al., 2004; Ootsubo et al., 
2007a). 

Aside from the crystals, Harker et al. (2007) find that 
the Deep Impact ejecta was rich in amorphous pyroxene, in 
contrast with the amorphous olivine and amorphous car- 
bon dominated pre-impact coma in July 2005, but simi- 
lar to the higher dust production rate coma in May 2005. 
In the impact-induced coma, crystalline olivine is corre- 
lated with the amorphous pyroxene dust, suggesting these 
minerals originated from the same reservoir (the Deep Im- 
pact crater). Harker et al. (2007) also found a population 
of amorphous carbon grains traveling at a projected speed 
of 700 m s _1 , equivalent to the highest speed ejecta seen 
by other observers (see, e.g., Keller et al., 2007). Schultz et 
al. (2007) interpret high speed ejecta as originating from 
the upper 5 m of the nucleus, whereas slower ejecta is from 
deeper layers. Harker et al. (2007) and Kadono et al. (2007) 
conclude that the surface of the nucleus is rich in sub-/xm 
amorphous carbon grains. 

With the broad wavelength coverage of Spitzer (5- 
40 /xm), Lisse et al. (2006) were able to identify a number 
of minerals that have not been previously reported in 
cometary spectra, namely Fe-rich silicates, metal sulfides, 
phyllosilicates, and carbonates. To emphasize the benefit 
of the long wavelength coverage, we assembled a bar chart 
consisting of 9P/Tcmpel's dust components (Fig. 4). The 
figure highlights those minerals that in principle arc diffi- 



cult to constrain without the 5-40 ^m Spitzer spectra — 
these components add up to one-third of the total dust 
surface area. Lisse et al. (2006) find that 8% of all the 
silicates by emitting surface area are in the form of the 
Fe-rich hydrated mineral nontronite, and 5% of the total 
surface area is in the Mg-rich carbonate magnesite, or the 
Fe-rich carbonate siderite. In §5.1, we discuss implications 
of the presence of hydrated minerals on our knowledge of 
comet nuclei. 

Harker et al. (2007) and Ootsubo et al. (2007b) assumed 
a Mg-rich composition for their crystalline silicate grains, 
which allowed adequate fits to most of their spectra (\ 2 ~ 
1.0). Lisse et al. (2006) simultaneously fit both Mg-rich and 
Fe-rich silicate minerals to their spectra and find significant 
amounts of Fe-rich silicate crystals. By number, the olivine 
Fe-/Mg-rich ratio is 0.3, and the pyroxene Fe-/Mg-rich ra- 
tio is 3.1. The difference between the Fe-/Mg-rich ratio in 
these two major silicate groups is intriguing. Both comets 
9P and Hale-Bopp have high crystalline silicate fractions, 
but their specific minerals are different. Only the Mg-rich 
species were present in comet Hale-Bopp (Wooden et al., 
1999; Bradley et al., 1999; Harker et al., 2002). Further- 
more, an 11.25 /im peak in the spectra of other OC comets 
indicates Mg-rich grains (Hanner et al., 1994). As discussed 
below (§5.2), if aqueous water was not present in the comet 
nucleus (§5.1), then the Fe content of silicate minerals may 
be a gauge of the amount of water vapor present in the 
solar nebula during silicate crystal formation. Future ob- 
servations and analysis of crystals in Jupiter-family and 
other comets from Spitzer and SOFIA are vital to our un- 
derstanding of the origin and evolution of dust in our solar 
system. 

5. The Implications of Hydrated Minerals and 
Fe-Rich Crystals on Dust Formation and Comet 
Evolution 

5.1. Hydrated Minerals 

The suggestion that Fe-rich crystalline silicates, phyl- 
losilicates, and carbonates are components of comet dust 
has consequences on our understanding of dust formation 
and cometary evolution. Significant amounts of liquid wa- 
ter would be required if these minerals were to have formed 
in 9P/Tcmpcl's nucleus (Lisse et al., 2006). This would up- 
set the common understanding that the interiors of comets 
remained at low temperatures after formation (Meech and 
Svoreh, 2004). Carbonates are thought to form by aqueous 
alteration inside asteroids (Armstrong et al., 1982; Zolen- 
sky and Browning, 1997), because under steady-state con- 
ditions in the PPD mid-plane, the nebular concentration 
of C0 2 is too low (Lewis et al., 1979). The (anhydrous) 
mineral olivine could have been aqueously altered to be- 
come the phyllosilicate serpentine in water-rich shocks in 
the solar nebula phase Ciesla et al. (2003) ; Fegley and Prinn 
(1989). In contrast to the relatively simple incorporation of 
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the OH— radical into olivine to form serpentine, the forma- 
tion of smectite requires mobilization of refractory element 
cations including Ca, Mg, and Si. The gas transport of re- 
fractory elements requires very high temperatures and this 
seems inconsistent with a shock-wave producing low tem- 
perature mineral phases such as phyllosilicates (A. Krot, 
private communication). Hence, phyllosilicates other than 
serpentine probably formed in the interiors of asteroids 
where pressures were higher than in the nebula (Alexander 
et al., 1989; Brearley, 2006; Scott and Krot, 2005). 

To aqueously alter minerals in the comet nucleus, some 
mechanism must warm the comet interior enough to melt 
or sublime water ice. Deeper than a few km, heating by 26 A1 
is sufficient to crystallize amorphous water ice, but much 
of the excess heat is quickly dissipated due the the higher 
heat conductivity of crystalline ice (about 20 times the heat 
conductivity of amorphous ice) . Only in comets with radii 
greater than 10 km can the temperature continue to rise to 
allow for liquid water to form (Prialnik and Podolak, 1995, 
1999; Prialnik et al., 2008; Merk and Prialnik, 2006; McK- 
innon et al., 2008). The interior may be high enough to 
form phyllosilicates by aqueous alteration of silicate min- 
eral grains in comet dust, but not hot enough to dry out 
the phyllosilicates to form FeO-rich crystalline silicates, as 
has been suggested to occur in asteroids (Krot et al., 1995, 
2000). In addition to phyllosilicates, other minerals that 
characterize parent-body aqueous alteration include hy- 
droxides, hydratcd sulfides, sulfates, oxides, and carbonates 
(Zolcnsky and Browning, 1997). These secondary minerals 
are notably absent from comet lP/Hallcy (Jessberger et al., 
1988; Jessberger, 1999) and Stardust (Brownlee et al., 2006; 
Zolensky et al., 2006, 2007). The Lisse et al. (2006) analysis 
of comet 9P suggests low concentrations of carbonates and 
the phyllosilicatc nontronite (a smectite), at the «5-10% 
level in comet 9P, but the identification of these minerals in 
spectra of other comets is still under study (Woodward et 
al., 2007; Bockelee-Morvan et al., 2007; Lisse et al., 2007). 
In an analysis of spectra of comet Hale-Bopp, (Wooden et 
al., 1999) set a tight constraint of < 1% on the phyllosilicate 
montmorilonite (a smectite) . Note that radiogenic temper- 
atures in comet nuclei are not sufficiently high, nor the con- 
ditions "dry" enough, to anneal and crystallize amorphous 
refractory silicate minerals in comet dust. As described be- 
fore, crystalline grains appear to have formed prior to the 
accretion of comet nuclei. 

5.2. Fe-Rich Silicate Crystals 

Since the discovery of crystalline silicates in comets 
(Campins and Ryan, 1989), it has been recognized that 
some comets exhibit mid-IR spectral features attributable 
to crystalline silicates, typically crystalline olivine. The 
identification of Mg-rich crystalline olivine at 11.2 ^m is se- 
cured by the ISO SWS spectrum of comet Hale-Bopp that 
shows the distinct longer wavelength features (Crovisier 
et al., 2000). On the other hand, some OC comets do not 



exhibit crystalline silicate emission features, although how 
much of this dichotomy is due to observational uncertain- 
ties will not be known until an analysis of a large number 
of comets is presented. Crystalline silicates are found in 
such high abundances in some comets that the crystals 
cannot be relics from the interstellar medium, where the 
crystalline fraction is < 2.2% (Kemper et al., 2004, 2005). 
For example, the crystalline-to-amorphous sub-/im mass 
ratio in the coma of the OC comet Hale-Bopp was 1.5- 
3.7 (Harker et al., 2002, 2004a). Recent observations have 
shown that JFCs also have crystals in significant frac- 
tional abundances: —0.3 for Deep Impact-induced inner 
coma (3") of 9P/Tcmpcl (Harker et al., 2005, 2007; Sugita 
et al., 2005) and —0.8 for the Deep-Impact larger coma 
(10" Lisse et al., 2006), -0.9 for 78P/Gchrels (Watanabe 
et al., 2005), and qualitatively described as less than in 
Hale-Bopp is 73P-B/Schwassmann- Wachmann and 73P- 
C/Schwassmann- Wachmann (Harker et al., 2006b; Sitko 
et al., 2006d). 

Exactly what conditions led to the formation of crys- 
talline silicates and their transport to the comet forming 
zone is not known, but the high Mg-content of crystalline 
silicates found in comets requires significant processing of 
proto-solar nebula amorphous silicates (see a review by 
Wooden et al. , 2007) . Based on a legacy of solar nebula con- 
densation models developed to explain minerals in chon- 
drules and chondrites, Mg-Fe amorphous silicates evapo- 
rated and re-condensed at — 1450 K as extremely Mg-rich 
silicate crystals and separate Fe grains in what is thought 
to be a typical "dry" solar nebula (Krot et al., 2000; Wein- 
bruch et al., 2000). In addition, annealing, the process of 
heating an amorphous grain enough to crystallize it, may 
have been an important process. Based on a wide variety of 
experiments designed to create amorphous silicates and an- 
neal them into crystals at —1000 K (see a review by Wooden 
et al., 2005), Mg-rich amorphous silicates anneal to Mg-rich 
crystalline silicates, Fe-rich amorphous silicates anneal to 
Fe-rich crystalline silicates, and Mg-Fe amorphous silicates 
anneal to Fe-bearing crystalline silicates (Brownlee et al., 
2005) but the Fe-content in the crystals begins low (0.1) and 
increases (to 0.2) as the crystallization completes (Murata 
et al., 2007). The annealing products and activation ener- 
gies somewhat depend on how the amorphous silicate start- 
ing materials are made in the laboratory (see the review by 
Wooden et al., 2005; Murata et al., 2007). The major issue 
with annealing to create cometary crystals is that Mg-rich 
amorphous silicates are not seen in cometary interplane- 
tary dust particles or chondritic materials, so there is no 
reservoir with which to anneal to Mg-rich crystals. Instead, 
annealing of Mg-Fe amorphous grains to form Mg-rich crys- 
tals must occur in a low oxygen fugacity gas (i.e., when 
water is depleted relative to molecular hydrogen) where Fe 
can be reduced (removed by interdiffusion) (Davoisne et 
al., 2006; Wooden et al., 2007). 

The small number abundance of Fe-rich crystals in 
cometary materials can be explained either by annealing 
without Fe reduction (e.g., Brownlee et al., 2005) or by 
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late-time condensation from the gas phase in a wet nebula 
(H2O/H2 wlOO-1000 times enhanced over a 'dry' nebula; 
Palme and Fegley, 1990). Interpretation of Stardust sam- 
ples of comet 81P/Wild, where Fe-rich crystalline silicates 
are found in small numbers (Zolcnsky et al., 2007) but 
where secondary aqueous alteration products (e.g., layer 
lattice silicates) are not found (Zolcnsky ct al., 2006), 
the gas phase condensation of moderately Fe-enriched to 
Fc-rich silicates is a more consistent scenario than parent 
body aqueous alteration. Moreover, we suggest that the 
condensation of Fe-rich crystalline silicates from a hot, 
water-rich nebula at times late enough for 90% of the Mg 
to have condensed into solids (Palme and Fegley, 1990) 
may explain the presence of Fe-rich crystalline silicates in 
both comets and in some chondrites that otherwise do not 
show extensive signs of parent body aqueous alteration. 

As discussed above, Fe-rich crystalline silicates can 
condense from the gas phase when the oxygen fugacity 
is high. In order to have both Mg-rich and Fe-rich crys- 
talline condensates, the water content must be variable. 
Water variability could come from inward migration of icy 
planctcsimals into the hot inner-zones of the PPD (Cuzzi 
and Zahnle, 2004). Therefore, the differences between the 
comets with Mg-rich crystals and comet 9P/Tempel with 
Fe-rich crystals may be a matter of the epoch of each 
comet's formation, and the efficiency of radial transport in 
and out of the comet forming zones. Studying the abun- 
dance of Fe-rich crystalline silicates in bodies that have 
not undergone internal aqueous alteration (i.e., comets) 
may aid our understanding of the transport of water from 
beyond the frost line into the inner solar nebula, and the 
subsequent enrichment of the outer solar system with 
condensates that formed in a water rich environment. 

6. Summary 

We reviewed ground-based and space-based mid-IR spec- 
tra of Jupiter-family comets taken over the past 25 years. 
These observations, including spectroscopy of 9P/Tempcl 
during the Deep Impact encounter and early Spitzer results, 
allow us to draw some conclusions on the nature of dust in 
JFCs: that, 1) JFCs have weak silicate emission features, 
roughly 20-25% over the continuum, similar to the lowest 
values observed in Oort Cloud comets; 2) the weak sili- 
cate emission could be due to low silicate content relative 
to other grain species (i.e., amorphous carbon and FeS), 
or due to the preponderance of large, compact grains (> 
1 /^m); 3) three JFCs show evidence for crystalline silicates 
in their mid-IR spectra (four if we include the Centaur- like 
29P/Schwassmann-Wachmann) — the existence of crystals, 
also found in OC comet spectra, indicate high temperature 
processes altered some dust in the proto-planetary disk; 4) 
JFC nuclei can have heterogeneous dust properties, the ev- 
idence for which is presented in studies of the Deep Impact 
encounter with comet 9P/Tcmpel; 5) the evidence for hy- 
drated minerals in the coma of 9P/Tempel is perplexing 



when we consider that comet interiors (for R < 10 km) 
likely remained well below the melting point of water ice; 
and, 6) the crystalline silicate forming zones of the solar 
nebula may have had a variable water content (i.e., vari- 
able oxygen fugacity) , which could lead to the apparent 
dichotomy of comets: comets that predominantly consist of 
Mg-rich crystals (e.g., comet Hale-Bopp), and comets with 
significant amounts of Fe-rich crystals (comet 9P/Tempel). 

We identify several opportunities with which we may im- 
prove our understanding of dust in JFCs: 1) increase the 
number of observed JFC mid-IR spectra; 2) continue to 
model the thermal emission spectra of JFCs in order to as- 
certain the nature of their weak silicate features; 3) con- 
tinue to compare models of light scattered by comet dust 
with thermal emission models to significantly constrain the 
structure of dust grains; and 4) employ the remarkable ca- 
pabilities of 8-10 m class telescopes and SOFIA to improve 
our understanding of the amount and mineralogy of silicate 
crystals in comets, especially their Mg and Fe contents. We 
hope to answer: Why do JFCs and some OC comets have 
weak silicate emission features? How common are Fe-rich 
silicates and hydratcd minerals in JFCs? Is the typical JFC 
heterogeneous in dust composition, and does that hetero- 
geneity arise from the primordial aggregation of planetesi- 
mals? 
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Wavelength (jam) 

Fig. 1. Model spectra of the grain components detect in post-Deep Impact observations of 9P/Tempel from Harker et al. (2007). The 
Gcmini-N/Michelle spectrum (black circles) is centered on the nucleus and is taken 74 min after impact. The thermal contribution of the 
nucleus is shown as an orange line. Superposed on the spectra is the total model spectral energy distribution (red line). The model is 
decomposed into its constituent dust components: amorphous pyroxene (blue line), amorphous olivine (cyan line), and crystalline olivine 
(green line). Notice the contrast between the broad spectral features of the amorphous silicates with the narrow resonances of crystalline 
olivine. 
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Fig. 2. The variation of selected resonance wavelength positions with Mg content in the 10 p.m region for crystalline olivine (Koike et al., 
2003) and crystalline ortho-pyroxene (Chihara ct al., 2001). Top — A post-Deep Impact Spitzer /IRS spectrum of comet 9P/Tempel (1 hr 
after impact) normalized by a scaled Planck function fit to the spectrum at 8 and 12 /im (as described in §2). The wavelength positions 
for silicate crystalline peaks observed in comets C/1995 Ol (Hale-Bopp) (Wooden et al., 1999), C/2001 Q4 (NEAT) (Wooden et al., 2004), 
and 9P/Tcmpcl (Harkcr ct al., 2007) spectra are marked with vertical dotted-lines. Center and Bottom — Filled symbols mark the central 
wavelengths of strong silicate resonances as they vary with Mg content, open symbols mark weaker silicate resonances. The y-axis values 
range from 0.0 (Mg-poor, Fe-rich) to 1.0 (Mg-rich, Fe-poor) corresponding to the x and y parameters in each mineral's chemical composition: 
[Mgj / ,Fci_ !/ ]2Si04 (olivine) and [Mg x ,Fei_2,]Si03 (pyroxene). 
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Fig. 3. Pre-impact spectrum of comet 9P/Tempel taken 03 July 2005 UT (gray line), and two spectra taken 0.64 (thick black line) and 20.5 hr 
(thin black line) after impact. The contribution from a 3.3 km nucleus has been removed. Notice the broad amorphous silicate emission at 
8-12 and 15-20 /im in the pre-impact spectra. Also note the strong crystalline resonances at 9 and 11 fim in the first post-impact spectrum. 
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Fig. 4. Grain components identified in the mid-IR in post-Deep Impact Spitzer spectra of comet 9P/Tcmpel (Lisse et al., 2006). In order 
to emphasize the importance of space-based infrared spectra, we have highlighted those components that have only been identified through 
reduced \ 2 nts to the Spitzer /IRS 5-35 /im observations. 
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Tabic 1: Mid-infrared spectroscopy of Jupiter-family comets. 



Comet 


Investigators 


Telescope Instrument 




1NOLCS 










Silicate Band 












Strength 




2P/Enckc 


Kelley et al. 2006 


Spitzer 


IRS 


N 


Silicates arc < 31% of the sub-/im dust 
mass (3ct). 


4P/Faye 


Hanner et al. 1996 


IRTF 


BASS 


1.25 




9P/Tempel 


Lynch et al. 1995 


IRAS 


LRS 


< 1.40 b 






Lisse et al. 2006 


Spitzer 


IRS 


Y 


We find the silicate band strength to 
be pre- .Deep Impact: 1.19-1.25 (Fig. 3); 
0.64 hr post-Deep Impact: 1.95-2.21 




Harker et al. 2007 


Gcmini-N 


Michelle 


Y 


We find the silicate band strength to be 
nre-Dera Imnact- Mav 2005 1 50 + 18 
July 2005 1.2 ± 0.2; 74 min post-Deep 
Impact: 2.02-2.16 




Ootsubo et al. 2007b 


Subaru 


COMICS 


Y 




lOP/Tempel 


Lynch ct al. 1995 


IRAS 


LRS 


< 1.40 b 




19P/Borrelly 


Hanner ct al. 1996 


IRTF 


BASS 


1.11 b 


1.13-1.19 after nucleus subtraction. 




Woodward et al. 2002 


IRTF 


HIFOGS 


N 




24P / Schaumasse 


Hanner et al. 1996 


IRTF 


BASS 


< 1.10 b 


Bare nucleus? 


26P/Grigg-Skjcllcrup 


Hanner ct al. 1984 


UKIRT 


UCL Spec. 


N 




29P/Scliwassmann-Wachmann a 


Stansberry et al. 2004 


Spitzer 


IRS 


1.10 b 


1.19 after nucleus subtraction. 


67P/Churyumov-Gerasimenko 


Hanner et al. 1985 
Kelley et al. 2006 


IRTF 

Spitzer 


FOGS 
IRS 


N 




69P/Taylor 


Sitko et al. 2004 


IRTF 


BASS 


1.23 b 


1.24-1.28 after nucleus subtraction. 


73P / Scliwassmann-Wachmann 


Harker et al. 2006a 


Gcmini-N 


Michelle 


1.15-1.25 


Fragments B and C 




Sitko et al. 2006a,b 


IRTF 


BASS 


1.18-1.25 


Fragments B and C 




Sitko ct al. 2006d 


Spitzer 


IRS 


1.33 


Fragment B 


78P/Gchrels 


Watanabc et al. 2005 


Subaru 


COMICS 


Y 




103P/Hartlcy 


Lynch ct al. 1998 


IRTF 


BASS 


1.15 b 


1.22 after nucleus subtraction. 




Colangcli et al. 1999 


ISO 


ISOPHOT 


N 


The uncertainties may be consistent 
with < 1.2 




Crovisier et al. 2000 


ISO 


ISOCAM 


1.20 b 


The nucleus is < 1% of the flux. 



a The Tisserand parameter of 29P is Tj = 2.98, which places the comet in the Jupiter-family class (Levison, 1996), yet in long timescalc 
integrations, the comet always remains between the orbits of Jupiter and Saturn (see, e.g., Carusi ct al., 1985) — a behavior that is similar 
to Centaurs. Since Centaurs and Jupiter-family comets have the same origin in today's solar system (the Kuiper Belt and Scattered Disk) 
their dust properties are related, and a comparison between 29P and JFCs is valid. 

b These silicate band strengths do not take into consideration the flux from the nucleus. For most entries, we list a corrected value in the 
notes column. 
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Table 2 

Coma (Fc) and nucleus (fjv) fluxes at 10 fim for our "median Jupiter-family comet" model. a 



Phase Fc 






A angle F N (6 = 0.20"; 


) b 0.8F N /(0.8F N + F c ) c 


(AU) (°) (mjy) (mjy) 


6 = 0.20" = 0.60" = 1.85" 


Afp = 110 cm, r h = 


= 1.8 AU, Rn = 1.8 


km 


0.8 158 9 


0.93 0.82 


0.60 


1.0 25 91 7 


0.91 0.77 


0.52 


1.2 31 60 6 


0.89 0.73 


0.46 


A/p = 430 cm d , r h = 


= 1.0 AU, Rn = 1.8 


; km 


0.3 81 1602 481 


0.73 0.47 


0.22 


0.6 72 488 240 


0.62 0.35 


0.15 


0.9 63 259 160 


0.56 0.30 


0.12 



a See §3.3 for a discussion of the model parameters. Note that Fc, Fff, and their ratio vary with wavelength. 
b The coma flux scales linearly with aperture angular radius, 9, for constant isotropic dust ejection (Eq. 2). 

c The aperture angular radius for the R 100 mode of three instruments: = 0.20" corresponds to Gemini/Michelle, 9 = 0.60" 

IRTF/MIRSI, and 8 = 1.85" to Spitzer/YRS. For simplicity, we assume a slit throughput of 80% for point sources (the nucleus). 

d To estimate the Afp at rh = 1.0 AU, we scaled the Afp at rh = 1.8 AU by a factor of 1.8 2 ' 3 as suggested by A'Hearn ct al. (1995). 



21 



